Diffusion between evolving interfaces 
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Diffusion in an evolving environment is studied by continuos-time Monte Carlo simulations. Diffu- 
sion is modelled by continuos-time random walkers on a lattice, in a dynamic environment provided 
by bubbles between two one-dimensional interfaces driven symmetrically towards each other. For 
one-dimensional random walkers constrained by the interfaces, the bubble size distribution domi- 
nates diffusion. For two-dimensional random walkers, it is also controlled by the topography and 
dynamics of the interfaces. The results of the one-dimensional case are recovered in the limit where 
the interfaces are strongly driven. Even with simple hard-core repulsion between the interfaces and 
the particles, diffusion is found to depend strongly on the details of the dynamical rules of particles 
close to the interfaces. Article reference: Journal of Physics: Condensed Matter 22, 465402 (2010). 

PACS numbers: 87.16. dp, 05.40.-a, 02.50.Cd, 66.10.cg, 64.60. De 



I. INTRODUCTION 

Diffusion phenomena are ubiquitous in nature, familiar 
examples ranging from heat conduction to osmosis. Of- 
ten diffusion occurs in a random or nonideal environment 
as it does, for example, in the presence of mobile or immo- 
bile (in the time scale of diffusion) lattice imperfections. 
Due to the complexity of real materials, transport has 
been considered within simplified theoretical frameworks, 
often utilizing the random-walk picture to describe some 
assumed underlying microscopy. Novel biophysical appli- 
cations can be expected to emerge for transport restricted 
by soft (and fluctuating) interfaces forming narrow or 
even nanoscale channels [l[. In particular the crossover 
from bulk-dominated to boundary-dominated diffusion is 
of considerable theoretical and experimental interest. 

The study of random walks in a random environment 
(RWRE) has a long history and since the results [H, H| 
from the 70's, as reviewed in Ref. [||, a vast amount of in- 
formation has been accumulated. This randomness has 
been considered to manifest itself as non-homogeneous 
transition rates 0-01 ■ Spatially, (frozen) transition rates 
can sometimes be described by random walkers like in 
the Sinai model @. In general, theoretical studies have 
mostly been limited to models, where the environment, 
including the possible geometric constraints, is either sta- 
tionary or fast compared with the jump rate of the walk- 
ers. The mathematical problem of the random walk in 
an uncorrelated fluctuating environment has been consid- 
ered in Ref. |{|. The asymptotics of diffusion in contin- 
uum under a random forcing in the presence of damping 
were analyzed in Ref. [Io| and diffusion in restricted ge- 
ometries with homogeneous transition rates was consid- 
ered in Ref. [ll[ • Two-species zero-range process [HI, [l3| 
with suitably chosen transition rates leads to dynamics 
which can be considered as a diffusing particle in an 



'Electronic address: 
t Electronic address: 



jaiiric.k.juntunen@jyu.fi 



j uha. t . mcrikoskiOj yu.fi 



evolving environment [14| . However, in existing stud- 
ies the focus has not been in (undriven) diffusion. In 
Ref. fl5j Sane et al. considered a situation, where parti- 
cles are immersed in a background fluid inside a narrow 
channel. From the point of view of a single particle in 
the dilute limit, this could be interpreted as diffusion in 
a dynamic environment. In that work the focus was on 
the transition from single file dynamics to Fickian diffu- 
sion. The existing studies on particle dynamics in the 
presence of interfaces are mainly for particles immersed 
in a driven liquid (l6j . A related problem, the influence 
of geometry fluctuations on lateral diffusion in biological 
systems was studied very recently in Ref. [l7j . where it 
is noted that geometry fluctuations at a finite scale can 
affect diffusion at all scales. 

In this Article, we consider diffusion in a dynamic re- 
stricting environment, inside open evolving 'bubbles' be- 
tween two interfaces. For this, we combine two most 
simple models, the solid-on-solid (SOS) model of inter- 
faces and the continuous-time random walk on a lattice, 
models that are widely used and known to describe well 
interface fluctuations and particle diffusion. To be more 
specific, we study diffusion on a lattice in the environ- 
ment produced by the dynamics of the BCSOS2-modcl 
introduced in Ref. [l8| . containing two non-intersecting 
interfaces driven against each other. Thus, due to the 
interface dynamics, the actual transition rates of a dif- 
fusing particle become dependent on time and position, 
when the jumps of the particle are possible only inside the 
'bubbles' between the interfaces. Diffusion in the hydro- 
dynamic limit will then depend on the dynamics of the 
bubbles, e.g. trough their growth and merging. We shall 
concentrate on cases, where the particles do not affect 
the motion of the interfaces so that the dynamics of the 
two interfaces, e.g. the bubble-size distributions and their 
correlations, are in principle known [l8l ] and thus, in ad- 
dition to simulations, analytical arguments can be devel- 
oped for various limits, which is the particular strength of 
the model. We use our model as a testing ground for var- 
ious ideas describing different regimes of behavior of the 
diffusion coefficient. In addition, we study the possible 
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consequences of various choices of microscopic dynamics 
on the interaction of the interfaces and the diffusing par- 
ticle. To reduce the dimension of the parameter space, 
the models we combine are relatively simple. However, 
most of our results are expected not to be dependent on 
the details of the model, but are characteristic of systems, 
where the size distribution of the bubbles and their dy- 
namics (at the bubble scale) inside a material or at an 
interface between two materials become the rate-limiting 
factor for diffusion (at the hydrodynamic scale). 

This paper is organized as follows. In Sec.|TT]we define 
the evolving environment provided by the two interfaces 
and the dynamic rules of the walker (particle) in detail. 
After that, in Sec. lIIIl we describe the algorithms needed 
for efficient simulation of the combined dynamics and the 
sampling of the main quantities. An experimentally ori- 
ented reader could first skip Sees. ITlllIIII and proceed to 
Sec. IIII D| where we briefly describe the parameters and 
their physical interpretation. Our results for the com- 
bined interface and particle dynamics are presented in 
Secs. lIVIfVl A concluding discussion of our results is given 
in Sec. EH 
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FIG. 1: A snapshot of the two interfaces hi and hi and the 
possible jumps of a diffusing particle with their attempt rates 
a, /3, 7 and v. Also the corresponding difference h- = hi — 
hi is shown. The jumps of the particle are possible only 
within the bubbles (open spaces) between the interface. The 
jump rate v is for the one-dimensional motion (see dynamics 
4 in text) depending only on h-. To simplify the discussion, 
when possible, in the text the various rate parameters are 
collectively referred to simply as /j, in such a way that the 
total rate for the particle to attempt a jump is /i. 



II. MODELS 

A. Model for the interface dynamics 

The evolving environment for diffusion is produced by 
the dynamics of the BCSOS2 model that we introduced 
and discussed in Ref. [ll]. Below we give only a brief 
description of the BCSOS2 model so that the dynamical 
rules for diffusion become well defined for the present 
study. The BCSOS2 model is constructed by letting two 
one-dimensional BCSOS interfaces [2(J interact with each 
others. 

The location or the 'height' of a single BCSOS interface 
is described by a function hi(x, t) such that, for every site 
x = 1, L, 

hi{x + l,t) - hi(x,t) = ±1, (1) 

where we without loss of generality assume that the possi- 
ble values of hi are integers. With this restriction on local 
configurations of the interfaces, only two kinds of pro- 
cesses, adsorption (hi locally increases) and desorption 
(hi locally decreases), are available. In our continuous- 
time model (see Sec. IIIip . the parameters pi and qi give 
the transition rates of adsorption and desorption events, 
respectively, for transitions allowed by the condition ([1} . 
In what follows, time and various rates are measured in 
units where Pi + qi = 1 . 

In the BCSOS2 model there are two interfaces, hi(x, t) 
and h2(x, £), such that hi(x, i) is an even (odd) integer for 
odd (even) values of x, see Fig. [TJ The coupling between 
the interfaces is produced by demanding that they cannot 
intersect: 

h 1 (x,t)>h 2 (x,t) for all x,t. (2) 

We also impose the periodic boundary conditions 
hi(x,t) = hi(x + L,t) for i = 1,2. In the full BCSOS2 
model, there are then four parameters [(pi, qi), (p 2 , 92)] 
defining the transition rates for the interfaces h\ and h 2l 
respectively. To further limit the parameter space, we 
shall restrict the discussion to the symmetric case p\ = qi 
and qi = P2 (see Ref. [TH) so that the behavior of the 
BCSOS2 system is defined by one parameter, the driving 
parameter / defined as 

f=p 2 /q 2 -l. (3) 

For large / the interfaces are strongly driven against each 
other and for / — > they become free. We also define the 
sum and difference processes defined via h± — h\ ± h 2 , 
where the sum process h + describes the wandering of the 
interfaces together and the difference process /i_ is posi- 
tive inside the bubbles and zero elsewhere [l8| . The non- 
crossing condition of Eq. is equivalent to h (x, t) > 0. 
The interfaces h± are then of the RSOS type [201, obey- 
ing h±(x+ 1, t) — h±(x, t) = —2, 0, +2. An example of hi 
and h 2 and the corresponding h^ configuration is shown 
in Fig. [D 
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B. Models for particle diffusion 

We consider a single point-sized particle diffusing be- 
tween the interfaces h\(x,t) and h,2(x,t) on a lattice 
[x,y). The lattice point coordinates in the horizontal 
direction are the same x = 1, ...,L as for the interface 
model above, again with periodic boundary conditions, 
c.f. Fig.Q] In the 'vertical' direction, the lattice is infinite 
and the coordinates are integers y = ...,—2,1,0,1,2,... 
and thus coincide with the possible values of hi and h?. 

We shall denote the location of the particle by (x p , y p ). 
The particle does not affect the dynamics of the interfaces 
but, if needed, a moving interface can push the particle 
the distance of one or two lattice units in the vertical 
direction such that the location of the particle also af- 
ter the change of the interface configuration satisfies the 
condition 



h 2 (x p ,t) <y p < hi(x p ,t). 



(4) 



These are moves of the particle forced by the interface 
motion. 

For diffusive moves of the particle, the following two 
rules are imposed in all cases: First, for a jump (x p ,y p ) — > 
(x' p , y' p ) to be possible, the product of the interface height 
differences on the departure site and the arrival site is 
non-zero: h-(x p ,t)h-(x' p ,t) > i.e. the channel for the 
jump between the interfaces must be open at both ends 
of the jump. Second, an attempted jump arriving outside 
the region bounded by the interfaces is blocked. In the 
actual dynamics, the direction of an attempted jump is 
chosen without any prior knowledge of the ability of the 
particle to perform the jump. 

For diffusion on the square lattice there are a few nat- 
ural choices for the possible particle jumps (x p ,y p ) — > 

(. x P i y p ) 

Dynamics m = 1: The most obvious case are the 
nearest-neighbor jumps such that the particle jumps in 
the horizontal direction (x p ,y p ) —> (x p ± l,y p ) with the 
attempt rate a and in the vertical direction (x pi y p ) — > 
(x p ,y p ± 1) with the attempt rate /3. We set a = (3 so 
that the total attempt rate of the particle is fx = Act. 

Dynamics m = 2: The particle jumps diagonally, i.e. 
(x p ,y p ) — > (x p ± l,y p ± 1) independently, with the at- 
tempt rate 7. The total attempt rate is then fj, = 47. 
This process is expected to be efficient on tilted sections 
of the interfaces like the rightmost part of the snapshot 
configuration in Fig. [TJ 

Dynamics m = 3: This is a combination of the jumps 
available in dynamics 1 and 2. In this work we chose 
a — ft = 7 so that /1 = 8a. 

In addition to the three models above, which we 
shall call two-dimensional particle dynamics, we con- 
sider simplified dynamical rules, which will be called one- 
dimensional: 

Dynamics m = 4: In this model only the x coordinate 
of the particle matters and the particle is allowed to per- 
form the jump (x p ,y p ) — > (x p ,y' p ) whenever the channel 
is open, i.e. h-(x p ,t)h-(x' p ,t) > 0. Then the effect of 



the interface dynamics on the possibility of the partic- 
ular jump is fully determined by h-(x,t) and diffusion 
is most directly controlled by the bubbles. In this case, 
in addition to the jump of the particle in the horizontal 
direction, there is, when needed, a move in the vertical 
direction over one lattice unit to keep the particle be- 
tween the interfaces such that Ii2(x p ,t) < y p < h\{x p ,t). 
The attempt rate of the jump in this case is denoted by 
v. The effect of forced moves on waiting times of the 
particle in the continuous-time dynamics is discussed in 
more detail in Sec. IIIIBI 



III. NUMERICAL METHODS 

A. Interface dynamics 

For the dynamics of the interfaces, which is the 
most time-consuming part of the numerics, we used the 



so-called iV-fold [21( algorithm in our continuous-time 
Monte Carlo simulations. In the iV-fold algorithm the 
possible transitions are divided into N classes according 
to their probabilities. After finding those classes, one 
finds all lattice points x, which belong to a certain class 
j. The next step is to calculate the set of time-dependent 
variables Qi = X)l=i n jk-Pjki where nj k is the number of 
those lattice points which belong to the class jk and Pj k 
is the probability associated with jk- The class j of the 
event, which will occur, is next determined by finding j 
such that Qj-i < R < Qj, where R is a random number 
with uniform distribution in the interval (0, Qn]- After 
finding the class, one randomly chooses a location (on h\ 
or /12) from this class. The waiting time for something 
to happen in the system consisting of the two interfaces 
then is Aij = — ln(i?i)/Qjv, where R\ is a random num- 
ber with uniform distribution in the interval [0, 1). 



B. Particle dynamics 

To restrict the dimension of the parameter space, in 
our simulations we only consider cases where all possible 
attempts (see Fig. [T]) of the particle to jump (within a 
given dynamics to) occur at the same rate as described 
in Sec. Ill Bl For simplicity and when possible, we shall 
denote by /i the total jump rate of a free particle, which is 
then determined by the parameters a, /3, 7 and v relevant 
in each case. In the presence of the interfaces this will 
become the total attempt rate such that with rate // the 
particle will try to jump in a direction chosen randomly 
from the list of allowed processes. The waiting time Atp 
after which the particle tries to jump is drawn form the 
exponential distribution Atp = — ln(i?2)//i, where R2 
is another uniformly distributed random number in the 
interval (0, 1]. 

In the combined model containing both the interfaces 
and the particle, the particle does not affect the dynamics 
of the interfaces, which will evolve as described in the 
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first paragraph above. To include the diffusing particle, 
we need two waiting times: the waiting time Atj of the 
interface system and the waiting time Atp of the particle. 
We then keep track of two times: First, the latest instant 
of real time i/, when there was a move in the interface 
system. Second, the latest instant of time tp, when there 
was either and attempt of the particle to jump or a forced 
move of the particle. Then, if ti + Ati < tp + Atp, the 
next event in the system will be a move of the interface, 
otherwise the next event will be an attempt to move the 
particle. 

After this, there still are in the continuous-time dy- 
namics two obvious ways to handle the forced moves (in 
dynamics m=l,2,3): (A) The clock of the particle re- 
mains intact in a forced move or (B) its waiting time 
Atp for the particle is updated after it. We shall con- 
sider both choices since they produce quite different re- 
sults and both can be physically justifiable from some 
microscopy. With this choice the dynamics of the com- 
bined model becomes defined. 



C. Sampling of main quantities 

The sampling of the observables was done with a con- 
stant time interval after reaching the steady state from 
an initial configuration consisting of two completely dis- 
ordered interfaces at a fixed distance from each other. 
With the A^-fold method reaching the steady state for 
interface configurations turns out not to be very difficult 
but especially for sampling of diffusion quite long runs 
were required. 

Typically of the order of 10 3 independent runs were 
performed for the interfaces in such a way that there were 
10 2 particles diffusing (independently of each other) be- 
tween the interfaces, the linear size of the simulation cell 
being L = 100 (for the finite-size scaling studies men- 
tioned in the text larger systems were used). The rea- 
son for this procedure is that the dynamics of the inter- 
faces even with the .TV-fold algorithm is computationally 
the most time-consuming part of the simulation, so we 
used each sequence of interface configurations to produce 
many particle trajectories, the total statistics thus being 
of the order of 10 5 . 

To characterize the statistical properties of the inter- 
faces hi and /12, we use their roughness or width (20| 
defined as 



\hi(x,t) — hi(x ± l,f)| = 1. In the notations for these 
quantities we suppress the dependence on system size L. 

The bubble size distribution per site was sampled at 
the same instant of times as the width of the interfaces. 
It gives the probability that a randomly chosen location 
x belongs to a bubble of size £, which is the length of 
the bubble in x direction. We shall be interested in 
bubble size distributions Pj(£,f), normalized such that 
only bubbles with j < £ < L are taken into account, 
J2f=j Pj(&> f) — 1j so that they are related by 



P j (£,f)=P (£,f) 



3-1 



t=0 



(6) 



In this £ = means that h-(x,i) = for site x. In 
the configurations, where the interfaces did not touch 
each other at all, the bubble size was recorded as a count 
in bin I = L. For certain purposes we also sample the 
bubble size distributions normalized for j < £ < L, which 
we denote by P*(£J), so that E^-P/ft/) = 1. In 
addition, we compute the (physical) average bubble size 
£(f) and also £*{}) defined as 



L-1 



L-1 



t(f) = J2 iP o^f)- (7) 



i=0 



e=o 



The mean waiting time for the change of the bubble size 
we shall denote by f(£,f). 

The diffusion coefficients, denoted here by I?obs, were 
determined in the long-time regime (\Ax(t)] 2 ) ; 



slope of the mean square displacement via [22| 
([Ax{t)} 2 ) ~2D ohs t, 



as a 



(8) 



where Ax(t) = x p (to +t) — x p (to) is the particle displace- 
ment during a time interval of length t. In the simu- 
lations, to reach the hydrodynamic regime not only for 
interfaces but for diffusion as well, we run the dynam- 
ics long enough so that the square root of ([Aa;(t)] 2 ) is 
much larger than the bubble size. Particle diffusion in the 
y direction, on the other hand, is in the long-time limit 
simply controlled by interface wandering and the effective 
sizc-dcpcndcnt diffusion coefficient related to that. 



W(f) = J(\h i (x,t)-h i (t))\< 



(5) 



Here i = 1,2 and hi{t) is the spatially averaged height 
of the interface configuration at time t, and the angle 
brackets denote ensemble average, i. e. average over inde- 
pendent simulations. The kink density k(f) is the density 
of those locations x, for which hi(x — l,t) ^ hi(x + l,t) 
for each of the interfaces i = 1, 2 separately. This def- 
inition follows from the fact that a BCSOS interface 
even in the flat state has an intrinsic roughness, because 
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D. Overview of the parameters and their meaning 

Our full model for diffusion restricted by fluctuating 
interfaces contais thus the following variable parameters: 

(i) The driving parameter /, defined in Eq. ^ of 
Sec. Ill A[ drives the interfaces towards each other. For 
interfaces to be coupled in the steady state we must have 
/ > 0. For increasing /, the channel for diffusion be- 
comes narrower and the lengths of the locally open paths, 
i.e. the bubble sizes, get smaller and the time scale of fluc- 
tuations gets longer. For / < the interfaces would be 
driven apart of each other and diffusion between them 
would become an unrestricted random walk. 

(ii) The second parameter is the generic jump at- 
tempt frequency /i of the diffusing particle. The possible 
jump directions (horizontal, vertical or diagonal, see also 
Fig. [T]) are controlled by the parameter m = 1,2,3,4 
described in Sec. Ill Bl The simplest choice is the one- 
dimensional jumps in the model with m = 4, which is 
a good starting point for analytical arguments for the 
scaling of diffusion with the model parameters. 

(iii) The particle jumps are restricted by the dynamic 
environment provided by the interfaces. The technical 
details of this are given in Sec. Ill Bl but there are no addi- 
tional parameters involved. However, there remains two 
physically reasonable ways to realize the occasions, when 
the interfaces would possibly move the particle. This is 
described in Sec. lIIIBl In scheme (A) the particle is con- 
sidered in such cases to move together with an interface 
with its diffusive 'clock' left intact and in (B) its 'clock' is 
updated after the interface move. The physical meaning 
of these choices is considered in Sec. IVII 



IV. RESULTS FOR INTERFACE DYNAMICS 



To obtain more detailed information, we computed T^ lb 
for j = 0,1,2 because the bubbles of size I = 0,1,2 
control diffusion for large /. 

In Fig. [2] we present the timescales T int , Ti at , T^ h 
describing the interface dynamics. The time T- lnt is de- 
termined by the interface configuration and its behavior 
is not monotonic as can be seen from Fig. [3J where a 
shallow minimum is observed. On the other hand, the 
roughness W as a function of /, shown in the same fig- 
ure, displays a dip at f w w 0.15. The non-monotonic 
behavior is explained by the entropic effects through the 
reduced configuration space available for the interfaces 
fl8| . The finite-size scaling of the dip position, f w (L), 
was studied in Ref. [IH, we only mention here the result 
f w ~ L" 1 / 3 for large L. In passing we note that such a 
deroughening due to interactions between interfaces has 
been experimentally observed in another context, see the 
articles in Ref. [19[ . It is evident from Fig. [5] that the 
point / = f w is also the crossing point for the dynamical 
properties of the interfaces. For / > f w , T mt increases 
rapidly and displays the asymptotic behavior Tj nt ~ / 
for large /, because the rate-limiting factor is the time 
required to create new bubbles, which is proportional to 
l/q-2 — (1 - <?2)/<72 - 1 = / for q 2 < in this limit. 

In Fig. [5] we also observe that T\ at ~ 1/f for / <C f w - 
This is controlled by the timescale of (finite) interfaces 
wandering apart from each other and then back together 
(effectively biased diffusion), which is proportional to 
1/(P2 - <?2) = 1/(2^2 - 1) ~ 1/f- For large / we have 
Ti a t 1/2, because the smallest bubbles I = 1 in this 
limit disappear with the rate 2qi —> 2. For increasing 
/, the time scale tS~1 increases without limit with the 
waiting time for a bubble to appear as T^ h ~ /, but 
Tu*u — > 1/2 by the same argument as for Ti at . Also the 



In this Section we first complement the study of 
Ref. [lH to characterize the time evolution of the coupled 
interfaces themselves (without a diffusing particle) in the 
stationary state. In the observed dynamics, a few time 
scales of interest can be monitored to gain insight also 
into the behavior of diffusion between the interfaces to be 
discussed in the next Section. The first time scale is the 
average time T; nt elapsed between consecutive changes in 
the interface configuration. The second one, denoted by 
Ti a t, is the average time scale over which a single loca- 
tion in x direction stays within a bubble i > in the x 
direction. Also the third one describes the behavior of 
the interface system only: It is the average waiting time 

(i) 

^bub f° r a change in the bubble size, 



L-l 



U) 



bub 



(9) 



where Pj (£, /) is the normalized bubble size distribution 
and f (£, /) is the corresponding average waiting time for 
something to happen for a bubble of size t, see Sec. MI Cl 




FIG. 2: Charasteristic timescales of the interfaces: the site- 
open time scale Ti at , the interface-change time scale Ti nt , the 
bubble waiting times T^ub f° r 3 ~ 0, 1,2 and the effective 
mobility time scale T mo b. 
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behavior of TJLt, has the same asymptotics, as can easily 
be seen by inspecting the one and only possible shape of 
bubbles of size I = 2, the rate- limiting factor being the 
shrinking of the bubble from its either end. 

The last time scale shown in Fig.[2]is the average wait- 
ing time for a change involving a bubble of size £ > 2 to 
occur in a given lattice site, which we denote by T mo b, 
since it is related to configuration changes that change 
the effective mobility of a particle by increasing or de- 
creasing its possible range of motion. This differs from 

(2) 

T} i in that, for example, for a particle sitting at site x, 
where the sites belongs to a bubble with I < 1, i.e. for a 
particle stuck in a locally closed configuration, in averag- 
ing T mo b we count the time for the particle to be mobile 
again, i. e. for the site to be in a bubble with £ > 2 again. 
For a bubble with £ > 2, on the other hand, each change 
of its size will result in a greater or smaller effective mo- 
bility for the diffusing particle. For / ^> f w we have 
Tmob ~ / 2 , since the rate-limiting process contributing 
to it in this limit is a two-step process, where a bubble 
of size two becomes created starting from a configuration 
where there are no bubbles. 

To characterize the properties of the interface configu- 
rations in more detail, we present in Fig.|4]the mean bub- 
ble size £ and the kink density k. The non-monotonicity 
of the differently normalized (see the caption) mean bub- 
ble size £* is a consequence of the fact that for finite 
L the maximum size of a bubble is limited. Below the 
roughness dip, for /</«,, the probability of completely 
open bubbles (£ > L), neglected in the computation of 
£, rapidly increases for decreasing /, while the number 
of such bubbles is essentially zero for / > f w . In the 
vicinity of the dip, for / ta f w , we observe I ~ /~ 4 / 3 , but 
for large / it tends to £ ~ f~ 1 , as expected. The kink 
density has a clear minimum slightly above f w , which 
is consistent with the reduction of configuration space 




FIG. 4: The average bubble size £ and kink density k as a 
function of /. We show the bubble size normalized in two 
ways: the real £ where bubble sizes < £ < L and £* where 
< £ < L are taken into account. 



due to the interaction between the interfaces, resulting 
in more hill tops and valley bottoms, which are not kink 
sites. Note also that k(f — > 0) = k(f — s- oo), because in 
the latter limit the interfaces are bunched together and 
their dynamics then reduces to that of a single isolated 
interface Hal. 



V. RESULTS FOR PARTICLE DIFFUSION 

For a continuous-time unbiased random walk the one- 
dimensional mean square displacement is of the form [22j 



([Ax(t)}*) = (N s )(ty Aa 



(10) 
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FIG. 3: The Ti nt and the scaled interfacial roughness W. Note 
that the vertical scale is different for these quantities. For this 
system size, the roughness dip is observed at / = f w « 0.15. 



where (N s )(t) denotes the average number of jumps in 
time t, and o\ x is the variance of the displacement of 
individual jumps. In what follows we shall use Eq. (|10[) 
to justify a theoretical model for diffusion between the 
interfaces for / > f w . For the different models of par- 
ticle dynamics in the presence of the interfaces, i.e. the 
models m = 1, 4 defined in Sec. Ill Bl we shall report 
our results as a function of the drive parameter / and 
the total jump rate fi as the dimcnsionless ratios 



(11) 



where ( f , n) is the observed diffusion coefficient ac- 
cording to Eq. © and D^l(n) is that corresponding 
to the same intrinsic jump rates without the interfaces. 
These ratios then give the effect of the interfaces on dif- 
fusion, while /i is the total attempt rate of jumps for the 
given model in units of the total attempt rate of local 
interface configuration changes. For a free random walk 
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with a\ x = 1, combining Eqs. (fT7|) and (fTU|) gives the 
diffusion coefficient for free motion in x direction, for ex- 
ample, as -Df, 1 ^, = a = (m/4 and D^ c = v = (m/2. Wc 
shall begin our discussion with the simplified dynamics 
m = 4, because for it the generic features of diffusion 
between the interfaces become more transparent. 

A. Diffusion for the one-dimensional particle 
dynamics (model 4) 

In the case of slow particles the interface configuration 
will change many times between the jump attempts of 
the particle. In this case, the success ratio of the jump 
attempts can be evaluated from the bubble size distribu- 
tion P *(£, /) and it is 

9(f) = J2PZ(l,f) l -=^+PZ(LJ), (12) 

1=2 

because only within bubbles with £ > 2 the particles 
can move and with probability (£ — l)/£ is an attempted 
jump possible since the particle cannot jump out of the 
bubble so that 2 of the 2£ attempts are blocked by the 
bubble edges. The mean-field prediction for the diffusion 
coefficient is then 

D^Hf,^=9(f)D^). (13) 

The simulation results for the model m = 4 together with 
this mean-field approximation are presented in Fig. [5] 
For fi <C 1/Tint, i.e. for slow diffusion, the curves D^f, fi) 
follow the mean field prediction of Eq. (|T3"|) . 

In the inset of Fig. [5] we show a data collapse by using 
for / the scaling factor 

a(n) = v/l + c^i (14) 




FIG. 5: Simulation results (plotting symbols) for Da{^,, /) and 
the mean-field prediction for the small fj, limit (solid line). In 
the inset we show the data collapse discussed in text. 



with c = 1. This is an interpolation of the large [i behav- 
ior cr(/i) ~ /i 1 / 2 for (j> 1, with the interface and parti- 
cle time scales well separated, and the small fi behavior 
<t(/x) ~ 1 for which is exactly the mean- field result 

discussed above. With this scaling, our simulation data 
for Di(f, /i) for / < f w is nicely collapsed onto the slow- 
particle curve defined by Eq. (|13|) and even better with 
the choice c w 1.1 (not shown), which we assign to finite- 
size effects. However, to obtain a good data collapse for 
f > fw we need a different scaling combination, 

"> ^/[l-cxp^r,)], (15) 
corresponding to an exponential clock with T2 = 

(2) 

r^ u ^(/ — > oo ) = 1/2. This scaling form results from the 
characteristic time scale of smallest bubbles allowing dif- 
fusion of particles in this regime. This is controlled by the 
driving parameter / such that for / 3> 1 there are mainly 
bubbles bubbles of size I = 1, but for diffusive jumps to 
take place bubbles of size £ = 2 are needed (with the 
time units chosen, once created, they stay open at least 
over one unit of time, cf. Sec.lIVp. This scaling is demon- 
strated in the loglog plot of Fig. [6] Because of the differ- 
ent rate-limiting mechanisms these scaling forms, below 
and above the dip region, are incompatible. 

Next wc consider the case of fast particles. If the par- 
ticle jump rate is fast compared to the interface dynam- 
ics, an adiabatic approximation can be done as follows. 
The particle is then trapped inside a bubble and its loca- 
tion becomes uniformly distributed in the timescale of the 
bubble dynamics and the (effective) location of the parti- 



ffi /i=0.001 




f 



FIG. 6: Data collapse for the large / behavior of D^ s by 
using the scaling form of Eq. (|15|) . The plotting symbols show 
the simulation results and the solid lines are the mean-field 
approximation for small of Eq. (|13|) with / scaled by the 
factor a(fj,) given in Eq. (|14[) . The dashed line is the adiabatic 
approximation of Eq. (|16[) and the dotted line is the finite- 
size-corrected infinite-rate approximation described in text. 
In the inset we show the comparison of the behavior of the 
jump- length factor B 2 and the mean squared bubble size I 2 . 
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cle can change only when the size of the bubble changes. 
The length of an effective particle jump via this mecha- 
nism is approximately the length of the displacement of 
the center-of-mass location b of the bubble. By applying 
Eqs. (|8I10[) we then obtain 

D a ^J) = \ 7 ^—B\f), (16) 

^ mob 

where T mo b is the time scale of bubble motion and B 2 (f) 
is the mean-square displacement (per jump) of the bubble 
obtained from 

B q (f) = (|6ncw - &old| 9 ), (17) 

with q = 2. Here 6 n ew and & id are the locations of a bub- 
ble before and after the change in the bubble size. This 
approximation is shown by the dashed line in Fig. [SJ For 
large / the expected behavior D ~ f~ 2 is recovered, since 
the bubble jump length and thus B 2 becomes a constant 
for isolated bubbles and T mo b ~ f 2 ■ For />/«,, the bub- 
ble motion is dominated by increasing or decreasing the 
bubble size by one lattice unit and, as expected, for fast 
particles the adiabatic approximation works well. For 
/ < fwj merging and dissociation of bubbles becomes 
important. We also tested the choice q = 1, with the 
mean jump length squared, giving less weight to long 
jumps and thus giving a smaller diffusion coefficient for 
small / (for large / the smallest jumps dominate in any 
case), but the approximation is not essentially better. 

The dotted curve just below the adiabatic approxima- 
tion shows an approximation obtained by assuming an in- 
finite jump rate of the particle such that after each inter- 
face configuration change a new location for the particle 
is drawn evenly distributed inside a bubble, but such that 
the particle displacement in the case of a completely open 
interface configuration is restricted by the system size. A 
possible way to extend this approximation for smaller / 
would be to utilize the known form of the bubble-size 
distribution for large I. This way we obtain a reduction 
of the diffusion coefficient that seems to work around 
/ = fw but apparently fails for small /, c.f. the behavior 
of £* in Fig. |4j Another natural approximation would be 
to consider the size of the bubble as defining an effec- 
tive mobility of the diffusing particle. Qualitatively, the 
behavior of B 2 and the mean squared bubble size £ 2 are 
quite similar, but with incompatible limits for small and 
large /. Unlike in Eq. (fT6"|) for B 2 , it turns out to be 
difficult to assign a natural rate factor to I 2 to develop a 
reasonable approximation for the diffusion coefficient. 

B. Diffusion for the two-dimensional particle 
dynamics (models 1,2,3) 

We next consider the effect of the 'microscopic' dynam- 
ics on diffusion. First, in Fig. [7] we show the diffusion co- 
efficient for dynamics m = 1 with the pushes of the parti- 
cle by the interfaces handled according to clock-updating 




f 



FIG. 7: Diffusion coefficient D\ for the particle clock updating 
scheme A. For comparison we show also by the full line the 
mean-field approximation, c.f. Fig. [5] In the inset we show 
the ratio D 2 (/x/2)/ D\(p). 
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FIG. 8: Diffusion coefficient Di for the particle clock updating 
scheme B. The full curves show the corresponding data for 
scheme A. 

scheme (A) and in Fig. [S] according to scheme (B), see 
Sec. Ill Bl for details. Now it is not just the length of the 
bubble but also the shape of it that affects diffusion. In 
(A) the particle waiting time is updated after each jump 
attempt, in (B) it is updated also after each push by the 
interfaces. The difference between (A) and (B) appears 
for long particle waiting times (for small /i) in the regime 
where there are frequent pushes, i.e. the drive / needs to 
be strong enough, but not so strong that diffusion would 
be completely dominated by the motion of small bubbles 
as discussed above. The coincidence of the waiting times 
of the interface motion and the particle jump attempts 
produces at a finite value of / in (B) behavior that looks 
like the suppression of diffusion sometimes observed close 
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to phase transitions. The dip in D\ appears close to the 
dip in the kink density k (see the inset of Fig. where 
in the local interface configurations there are more sites 
with possible interface configuration changes (at a kink 
site the interface is locally frozen) and thus more pushes. 
Other variants of the jump rate modification are con- 
ceivable, e.g. physically it might happen that the con- 
tact with the interfaces could boost the particle jump 
rates, but in the region of the parameter space, where 
bubble dynamics is the rate-limiting factor, it would not 
essentially change diffusion from what what is observed 
in scheme (A). 

The corresponding mean-field approximation of 
Eq. (fT3|) shown in Fig. [7] by the full curve is presented 
to help the comparison between m = 1 and m = 4. As 
evident from the comparison of Fig. [7] and Fig.[5j scaling 
by the factor a(fi) of Eq. [TJ] would not yield a very good 
data collapse for / <C f w in the case of two-dimensional 
particle dynamics. For particles fast both in the hori- 
zontal and vertical direction, the blocking of diffusion by 
'collisions' with the interfaces is quite efficient and leads 
to a considerable reduction of D\ in this regime. Note 
that even for small / we have D\ < 1 for fast particles, 
since particle diffusion in the vertical direction is faster 
than that of the interfaces for any / > 0. This in part 
explains the spreading of the curves -Di(/) for / < f w . 
In the inset of Fig. [7] we show also the difference between 
the particle dynamics m = 1 and m = 2 with the clock 
updating scheme A (concerning dynamics m = 3, we 
find that D\ < D$ < D2). For reasonable comparison, 
the jump rates for each dynamics have been chosen such 
that the total jump rates in the x direction match, thus 
we plot 1)2(^/2)/ Di(fi). Due to the limited accuracy of 
D m , the ratio D2/D1 becomes somewhat noisy, but a 
few general observations can be made. A difference first 
develops for f > f w , because diffusion becomes more re- 
stricted by the interfaces within bubbles of the type seen 
on the right (x = 17, 20) in the snapshot of Fig. [T]bc- 
coming prevalent. In such interface configurations, many 
of the 'horizontal' jumps (for m = 1) become blocked and 
diffusion along the narrow channel requires also 'vertical' 
jumps, while the diagonal jumps (for m = 2) arc more 
effective for particle transport. However, for / 3> f w dif- 
fusion again becomes dominated by bubble motion the 
way it was for m = 4, and D m ~ f~ 2 for m = 1,2,3 
for both particle clock updating schemes (A) and (B) as 
seen in Fig. [8] 



VI. DISCUSSION 

To summarize, we have considered one and two di- 
mensional continuous time random walkers constrained 
between two evolving interfaces symmetrically driven 
towards each other. A surprisingly complicated phe- 
nomenology appears. First, in the interface model itself 
there is a dip in the interface roughness at a finite value 
f w of the parameter / describing the drive of the inter- 



faces against each other [18[ , with different kinds of inter- 
face and bubble dynamics below and above the dip. The 
crossover from bubble-dominated (interface-dominated) 
to almost free diffusion is controlled by the relative jump 
rate of the particle and its interplay with the rate of the 
interface time evolution. In the analysis of this crossover, 
simple scaling arguments for the two regimes are incom- 
patible, both physically and formally, so that the full be- 
havior of diffusion cannot be described by a single scaling 
form. This can be expected to be a generic property of 
transport restricted by fluctuating interfaces. 

Furthermore, diffusion was found to depend on the mi- 
croscopic details of the interaction between the interfaces 
and the diffusing particles. In particular, its immediate 
effect on the waiting time of particle jumps was shown to 
be considerable, especially for dynamics physically cor- 
responding to diffusion on a lattice in the large-friction 
limit, which can be realized at domain boundaries in 
adsorption systems. In a spatial continuum this effect 
would be absent, as would be the roughness dip. How- 
ever, the finite-size case, with the underlying microscopic 
structure not washed out in coarse graining, can be of 
interest in its own right in nanoscalc applications. Then 
also the size of the channel for the particles will induce a 
relevant length scale, in addition to the length scale re- 
lated to the diffusive jumps (lattice) and the one related 
to the environment (bubbles). 

In this work we have considered diffusion that is ef- 
fectively one-dimensional even if the diffusive jumps and 
interactions between the tracer particles and the inter- 
faces result from two-dimensional dynamics. The envi- 
ronment for the diffusing particles is then described by 
a chain of bubbles, where as in higher dimensions more 
complicated topologies (networks) would arise, in some 
cases leading to a percolation problem. In such stud- 
ies, like in the present one, a considerable problem is the 
wide gap between the 'microscopic' timcscales related to 
the dynamics of the particles and the environment, and 
the 'hydrodynamic' timescale corresponding to diffusion 
over length scales larger than any structures in the en- 
vironment experienced by the particles. We hope our 
work would inspire further theoretical and experimental 
studies of diffusion in evolving and constraining environ- 
ments. 
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